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Abstract. We report simultaneous multi-frequency 
observations of the blazar PG 1553+113, that were 
carried out in March-April 2008. Optical, X-ray, 
high-energy (HE: > 100 Me V) 7-ray, and very-high- 
energy (VHE: > 100 GeV) 7-ray data were obtained 
with the KVA, REM, RossiXTE/ASM, AGILE and 
MAGIC telescopes. This is the first simultaneous 
broad-band (i.e., HE+VHE) 7-ray observation of 
a blazar. The source spectral energy distribution 
derived combining these data shows the usual double- 
peak shape, and is interpreted in the framework of 
a synchrotron-self-Compton model. 
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I. Introduction 

The transformation of gravitational energy from an 
accreation disk around a supermassive black hole into ra- 
diation is commonly believed to be the underlying cause 
of emission in Active Galactic Nuclei (AGNs). Further- 
more, much of the emission is relativistically beamed 
perpendicular to the disk by a mechanism that, although 
not completely understood yet, is most likely related 
to the focusing properties of the rotating, fully-ionized 
accretion disk (8). It is believed that the viewing angle 
of the observer determines the observed phenomenology 
of AGNs ll33l . The AGNs whose relativistic plasma jets 
point towards the observer are called blazars. The blazar 
class includes Flat-Spectrum Radio Quasars (FSRQs) 
and BLLac objects, the main difference between the 
two classes being their optical emission lines, which are 
strong and quasar-like for FSRQs and weak or absent in 
BL Lacs. 

The overall (radio through 7-rays) spectral energy 
distribution (SED) of blazars usually shows two broad 
non-thermal continuum peaks. For high-energy-peaked 



BL Lac objects (HBLs), the first peak is in the UV/soft- 
X-ray bands [as opposed to IR/optical for low-energy- 
peaked BLLac objects (LBLs)] whereas the second 
peak is in multi-GeV band (multi-MeV for LBLs). 
The low-energy peak is most commonly believed to be 
synchrotron radiation from a non-thermal population of 
electrons moving in a tangled magnetic field threading 
the emitting plasma region - the latter is modeled as a 
sperical blob that moves relativistically along the jet. The 
second peak forms by Compton upscattering of lower- 
energy photons by the same population of electrons that 
is responsible for the synchrotron emission. The lower- 
energy photons can originate either from the synchrotron 
radiation constituting the low-energy bump (Synchrotron 
Self Compton: SSC), or from outside the relativistic 
plasma blob (External Compton: an external source of 
'seed' photons could be the accretion disk [ 1 1 1 and/or 
broadline region [32]). 

Blazars often show violent flux variability, that may 
or may not appear correlated in the various energy 
bands. Simultaneous observations are then crucial to 
understand any patterns of temporal variability as well 
as the emission processes. 

The high-energy-peaked BL Lac source PG 1553+1 13 
was firmly detected at E > 200 GeV by MAGIC at an 
8.8cr significance level, based on the data from 2005 and 
2006 H. The source had also been tentatively detected 
at VHEs by H.E.S.S., at a 4cr significance level (5.3<r 
using low-energy threshold analysis: J2), which was 
confirmed later using 2005 and 2006 data O. The lack 
of detection of spectral lines (neither in emission nor 
in absorption) in the optical spectrum of PG 1553+113 
makes it impossible to measure its redshift directly 
lfl5l . However, an ESO-VLT spectroscopic survey of 
unknown-redshift BLLac objects suggests z > 0.09 
QUI . On the other hand, the absence of a break in the 
intrinsic VHE 7-ray spectrum may sugges z < 0.42 l27l . 
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II. VHE 7-RAYS : MAGIC OBSERVATIONS 

The MAGIC Telescope 0, (9) is a latest genera- 
tion Imaging Atmospheric Cherenkov Telescope (IACT) 
at La Palma, Canary Islands, Spain (28.3°N, 17.8°W, 
2240 m a.s.l.). Thanks to its low threshold of 50GeV 
1 5 1, MAGIC is well suited for multi-frequency obser- 
vations, together with the instruments operating in the 
GeV range. The parabolic shaped reflector, resulting 
in a total mirror area of 236 m 2 , allows MAGIC to 
collect Cherenkov light from particle showers initiated 
by gammas or other particles in the atmosphere and 
focus it onto a multi-pixel camera, composed of 577 
photomultipliers. The total field of view of the camera 
is 3.5°. The incident light pulses are converted into 
optical signals and transmitted via optical fiber to a 3- 
level trigger system. The selected events are digitized by 
2 Gsamples/second Flash ADCs ||2TI . With a statistical 
analysis of the recorded light distribution and the ori- 
entation of the shower image in the camera, the energy 
of the primary particle and its incoming direction are 
reconstructed. 

The MAGIC observations used in this paper were 
performed on 16-18 March and 13, 28-30 April 2008. 
The zenith angle of the data set ranges from 18 to 36 
degrees. Observations were performed in wobble mode 
ifTTl . where the object was observed at 0.4 degree offset 
from the camera center in opposite directions every 20 
minutes. After the quality selection of the data the total 
effective observation time is 7.18 hours. 

An automatic analysis pipeline (see [13]) was used 
to process the data. After an absolute calibration with 
muons |20l . and an absolute mispointing correction, 
noise subtraction and background reduction is achieved 
using the charge distribution and the arrival time in- 
formation of the pulses of the pixels [6|. Three OFF 
regions, defined in the same data set symmetrically 
to the ON region with respect to the camera center, 
were used to determine the remaining background. To 
select the gamma-like events, a dynamical cut in Area 
(Area=7r WIDTH LENGTH) versus SIZE and a cut in d 
were applied Q. More details of the cuts can be found in 

ED. 

For the calculation of the differential energy spectrum, 
a loose cut in Area was selected to ensure that more 
than 90% of the MC gammas after the image cleaning 
survive the cut. To study the dependence of the spectral 
shape on the cut efficiency, a different cut in Area with 
cut efficiencies of between 50% and 95% for the entire 
energy range was applied. The gray area in Fig.l is 
the result of this study. For data affected by calima 
(i.e. windblown sand dust from Sahara in an air layer 
between 1.5 km and 5.5 km a.s.l. causing absorption of 
the Cherenkov light), a correction has been applied. The 
method is described in detail in fl4l . 

'Width, Length and Size are the image parameters in [22] . 



III. HE 7-RAYS: AGILE OBSERVATIONS 

The Gamma-Ray Imaging Detector (GRID: 30MeV- 
30 GeV) on board the spaceborne HE-7-ray telescope 
AGILE @ observed PG 1553+113 in three different time 
periods: 16-21 March, 25-30 March and 10-30 April 
2008. The GRID data were analyzed using the AGILE 
standard pipeline (see [37] for a detailed description of 
the AGILE data reduction), with a bin size of 0.25° x 
0.25° for E > 100 MeV Only events flagged as con- 
firmed 7-rays and not recorded while the satellite crossed 
the South Atlantic Anomaly were accepted. We also 
rejected all the events with reconstructed direction within 
10° from the Earth limb, thus reducing contamination 
from Earth's 7-ray albedo. The total effective observation 
time, after the selection is 268 hours. PG 1553+113, 
observed at about 50 degrees off-axis with respect to 
the boresight, was not detected by the GRID at a 
significance level > 3 a and therefore an upper limit 
at 95% confidence level was calculated. Considering 
that AGILE has a higher particle background at very 
high off-axis angles, we calculated also the upper limit 
selecting only photons with E 200 MeV in order to 
minimize the possible contamination at low energies. 
The summary of the AGILE observations and the results 
of the analysis are reported in Table 1. During March- 
April 2008, PG 1553+113 was outside the field of view 
of SuperAGILE, the hard X-ray (20-60 keV) imager on 
board AGILE ifTBI . thus no SuperAGILE data on this 
source are available from this period. 

IV. X-RAYS: RXTE/ASM OBSERVATIONS 

The All Sky Monitor (ASM) on board the Rossi X-ray 
Timing Explorer (i?XTE) satellite consists of three wide 
angle scanning shadow cameras. The cameras, mounted 
on a rotating drive assemble can cover almost 70% of the 
sky in every 1.5 hours ]25'|. F° r this campaign, we used 
ASM data taken from 1 March through 31 May, 2008. 
The mean measured flux of PG 1553+113 is showed in 
Table 2. 

V. Optical data 

A. KVA observations 

The KVA (Kungliga Vetenskaps Akademien) tele- 
scope is located at the Roque de los Muchachos Ob- 
servatory on La Palma island, Spain and is operated 
by Tuorla Observatory. The telescope is composed of 
a 0.6 m f/15 Cassegrain devoted to polarimetry, and 
a 0.35 m f/11 SCT auxiliary telescope for multicolour 
photometry. This telescope has been successfully oper- 
ated remotely since autumn 2003. The KVA is used for 
optical support observations for MAGIC by making R- 
band photometric observations typically one measure- 
ment per night per source. Photometric measurements 
of PG 1553+113 were made in differential mode, i.e. 
by obtaining CCD images of the target and calibrated 

2 Astro-rivelatore Gamma a Immagini LEggero; (34), [35]. 
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comparison stars in the same field of view (|18|; |fl9ll ; 

ED). 

B. REM observations 

REM (Rapid Eye Mount: a fast slewing robo- 
tized infrared telescope, |10|) acquired photometry of 
PG 1553+1 13 on April 18 and 25, and May 2, 2008 with 
all the available filters (VRIJHK). The data reduction 
followed standard procedures [12]. The resulting flux of 
the three observations is reported in Table 2. The optical 
magnitudes have been calibrated against the standard 
Landolt field El of PG1323, while the near infrared 
magnitudes were calibrated against the 2MASS catalog. 
The optical data taken during the night of 18 April are 
affected by moonlight and have not been used. For the 
purpose of SED construction, all magnitudes have been 
dereddened using the dust IR maps PP . 

VI. Results 

Analysing the VHE data, an excess of 415 gamma- 
like events, over 1835 normalized background events 
was found, corresponds to a significance of 8 a. The 
differential VHE photon spectrum of PG 1553+113, re- 
sulting from MAGIC data averaged over all observation 
nights, is plotted in Fig.l (filled circles). It can be well 
described by a power law: 



dN / E \(-3.4±o.i) 



dE v V200GcVy (1) 

in units of lO^cm^s^TeV- 1 with x 2 /DOF= 1.36/3. 

During the nights of MAGIC observations no signif- 
icant flux variability was found. The differential energy 
spectrum for the March data alone can be well described 
by: 



E 



(-3.5±0.2) 



^ = ( 1 - 9±0 - 4 )-( 2 00GeV; 
whereas the energy spectrum for the April data can be 
well described by: 

dN E \(-3.3±0.2) 

V200GeV/ 



§ = < 2 .i ± o. 4 >.( ; 



(2) 



(3) 



(same units as above). Because of the interaction of VHE 
7-rays with the extragalactic background light (EBL: 
the integrated light from all stars that ever formed; it 
spans the IR-UV frequency range), which leads to e + e~ 
pair production and ensuing attenuation of the VHE 7- 
ray flux, we computed the deabsorbed (i.e., intrinsic) 
fluxes using a specific Tow star formation model' of the 
EBL ll23l . assuming a source redshift of z — 0.3. The 
resulting deabsorbed points are represented as empty 
squares in Fig.l. 

The HE data reduction results from AGILE are sum- 
marized in Table 1 . The upper limits obtained by AGILE 
are consistent with the average flux observed by Fermi- 
LAT for this source in its first three months of operation 
IT). The upper limit obtained in the second period (25- 
30 March) is of the same order of that obtained in the 
first period. 




200 300 400 
E [GeV] 

Fig. I. The MAGIC measured spectrum of PG 1553+113 (filled 
circles). The EBL-corrected points are shown as empty squares.) 

TABLE I 

The upper limit calculated by AGILE data in two 

DIFFERENT PERIODS. 



Time period 


Energy 


U.L. Flux [photons cm 2 s i ] 


16-21 March 


> lOOMeV 


56 X 10~ s 




> 200 MeV 


36 x 10~ 8 


10-30 April 


> lOOMeV 


34 x 10~ s 




> 200 MeV 


21 x 10~ 8 



The fluxes, and their corresponding effective photon 
frequencies, from the other telescopes used in this paper 
are reported in Table 2. 

TABLE II 

Effective frequencies of operation, and corresponding 
fluxes from pg 1553 + 1 13, from the various telescopes 
used in this paper. 



Instrument 


\og(v [Hz]) 


\og{v¥(y) [erg cm 2 s i ]) 


KVA 


14.63 


-10.2 


REM 


14.38 


-10.33 




14.27 


-10.34 




14.13 


-10.38 


RXTE/ASM 


18.03 


-10.3 



VII. Discussion 

The SED of PG 1553+113 is shown in Fig.2. The 
VHE and HE 7-ray points are from MAGIC and AGILE, 
repectively. The X-ray point, provided by RXTE/ASM, 
represents the average flux between March 1 and May 
31, 2008. The optical i?-band point, provided by the 
KVA telescope, is from March 18,19. 

In addition to these data, we also used an optical 
flux on April 18 from REM. To assess the soundness of 
the addition of this quasi-simultaneous data with other 
simultaneous data, we checked the optical variability of 
the source from April 13 to April 18 using KVA data, 
and found that the source was essentially stable, with a 
variability of less than 2% in flux. In order to see the 
variability in HE range from the upper limits of AGILE, 
we used flux points from the Fermi Gamma-ray Space 
Telescope (formerly GLAST; Abdo et al. 2009). 

We fit the resulting SED with a homogeneous one- 
zone SSC model |36| . In the frame of this model, the 
source is a spheric blob of plasma of radius R, moving 
with a Doppler factor <5 1 28 1 towards the observer at an 
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Log v [Hz] 

Fig. 2. The average SED of PG 1553+1 13. The filled triangles denote 
the REM data. The open circle represents the KVA data point. The 
open square is from HXTE/ASM. At HE 7-rays, the upper limit is 
from AGILE. The filled circles show the deabsorbed MAGIC data 
points. We also show the non-simultaneous flux from Fermi (butterfly 
symbol). 

angle 9 with respect to the line of sight and threaded 
with a tangled uniform magnetic field, B. The injected 
relativistic particle population is described as a broken 
power-law spectrum with normalization K, extending 
from 7 m ; n to 7 max and with indices n\ and n.2 below and 
above the break Lorentz factor 7b r . In the framework of 
this model, the SED can be described with the following 
set of parameters: 7 min = 1, Tbreak = 3 x 10 4 , 7 max = 
2xl0 5 , K = 0.5xl0 4 cm 3 , n x = 2, n 2 = 4, 5 = 0.7G, 
R = 1.3 X 10 16 cm, and S = 23. 

The Fermi and lowest-energy MAGIC data points 
do suggest some variability at HE 7-rays. Comparing 
the current SED with that based on earlier MAGIC 
observations @), some variability is clearly visible. This 
is due to a small flux variation in the X-ray range and a 
large flux variation in the VHE range. The optical flux 
does not show any significant variability. 

Our results suggest that the variablility of 
PG 1553+113 at different frequencies is highly 
time depended: hence, only a simultaneous multi-A 
monitoring over a large time span will give more 
information on the source. Relative to this fact, it 
is worth mentioning that the AGILE and MAGIC 
data presented here constitute the first simultaneous 
broad-band 7-ray observation (and ensuing SED) of 
any blazar. 
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